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Abstract. The anisotropic thermoelectric transport properties of bulk silicon 
strained in [lll]-direction were studied by detailed first-principles calculations 
focussing on a possible enhancement of the power factor. Electron as well as 
hole doping were examined in a broad doping and temperature range. At low 
temperature and low doping an enhancement of the power factor was obtained 
for compressive and tensile strain in the electron-doped case and for compressive 
strain in the hole-doped case. For the thermoelectrically more important high 
temperature and high doping regime a slight enhancement of the power factor 
was only found under small compressive strain with the power factor overall 
being robust against applied strain. To extend our findings the anisotropic 
thermoelectric transport of an [lll]-oricnted Si/Ge superlattice was investigated. 
Here, the cross-plane power factor under hole-doping was drastically suppressed 
due to quantum-well effects, while under electron-doping an enhanced power 
factor was found. With that, we state a figure of merit of ZT= 0.2 and ZT= 1.4 at 
T = 300 K and T = 900 K for the electron-doped [lll]-oriented Si/Ge superlattice. 
All results are discussed in terms of band structure features. 
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1. Introduction 

Thermoelectric phenomena were first described for metals by Seebeck at the beginning 
of the 19th century and revived by Ioffe in the late 1950's by the introduction of 
semiconductors to thermoelectric devices [1, 2]. However, since then thermoelectrics 
were restricted to a scientific and economic niche mainly due to their poor 
conversion efficiency [3, 4]. Nowadays emerging global need for energy production 
and conservation has intensified interest and research in more effective alternative 
energy technologies to reduce our dependence on fossil fuels. Contributing to this, 
thermoelectric devices could partially convert wasted heat into electricity by their 
ability to transform heat directly into electric current, and vice versa [5]. 

The thermoelectric conversion efficiency can be stated by the figure of merit 
(FOM) 

ZT = a f T, (1) 

where a is the electrical conductivity, S the thermopower, K e i and k p h are the electronic 
and phononic contribution to the thermal conductivity, respectively. The numerator 
of Eq. 1 is called power factor PF = aS 2 and characterizes the electric power output. 

While thermoelectric devices are extremely facile, have no moving parts, and 
do not produce greenhouse gases [6], two obstacles limit their applicability. The 
first, a low efficiency, could be challenged by nowadays nanostructurcd thermoelectrics 
enabling large values of ZT 1 [7, 8, 9]. As a second drawback, those materials are 
based on environmentally hazardous or rare lead, tellurium or selenium compounds 
and are therefore hard to integrate in semiconductor electronics. 

However, current research gained a tremendous progress in enabling silicon for 
thermoelectrics. Silicon is non-toxic, readily available, cheap and well integrated in 
present electronics infrastructure, so this might be a considerable leap forward. While 
silicon has been stated as inefficient thermoelectric in the past due to its enormous 
thermal conductivity [10], recent experimental and theoretical attempts revealed that 
nanostructuring could lead to thermoelectric efficiencies comparable to state of the 
art commercial thermoelectric materials [11, 12, 13, 14]. 

Besides the reduction of thermal conductivity, that is the denominator in Eq. 1, 
Koga et al showed in a seminal work [15, 16] that it should be possible to enhance the 
power factor, that is the numerator in Eq. 1. This concept of carrier pocket engineering 
uses the influence of strain to optimize the band structure of silicon and germanium 
based superlattices (SL) regarding their electronic transport. As a main result it was 
found, that the effect of the lattice strain at the Si/Ge interfaces is more relevant for 
strain in [111], than in [001]-direction, regarding a possible enhancement of the power 
factor. In fact, ZT = 1.25 and ZT = 0.98 at room temperature were predicted for 
strain-non-symmetrized and strain-symmetrized [lll]-oriented Si/Ge-SL, respectively, 
and the ZT values are shown to increase significantly at elevated temperatures [17]. 

While in a previous study we already concentrated on the influence of biaxial 
in-plane strain in [001]-direction on the thermoelectric properties of silicon [18], we 
will focus here on the influence of strain along the [111] -direction. For this purpose 
the paper is organized as follows. In section 2 we introduce our approach footing on 
first principle electronic structure calculations within density functional theory and 
transport calculations based on the solution of the linearized Boltzmann equation. By 
this means we start the discussion of the thermoelectric transport properties of bulk 
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silicon strained along the [111] -direction in dependence on strain and doping to gain 
inside into the physical mechanisms, which clearly differ from the [001]-strain case. 
The discussion is performed for electron- as well as for hole-doping in sec. 3.1 and 3.2, 
respectively. 

To extend the findings for the strained bulk silicon, in sec. 3.3 results for 
an exemplary Si/Gc-SL grown on Si in [111] -direction are presented. Here the 
influence of tensile strain in [11 Indirection, induced by the lattice mismatch at 
the Si/Ge interface, is investigated with respect to the thermoelectric transport in- 
plane and cross-plane the SL. Again, the temperature and doping dependence of 
the thermoelectric properties are discussed for electron- and hole-doping regarding 
a possible enhancement of the power factor. A further aspect will be the influence of 
structural relaxation and chemical composition on the transport properties. At the 
end of the paper, in sec. 3.4 insights into the FOM will be presented along with the 
electronic part of the thermal conductivity, to give a clue on optimal charge carrier 
concentrations to obtain the best FOM. 

While focussing our interest in the high-temperature thermoelectric application of 
strained silicon, our results in the room temperature regime could be of importance for 
the metal-oxide-semiconductor device community. Knowledge of the thermoelectric 
properties of silicon under strain could help to understand parasitic effects on the 
electronic transport in those structures. In this low temperature and low doping 
regime we confirm a remarkable influence of externally applied strain on the electrical 
transport under electron and hole doping. 

2. Methodology 

Our approach is based on two constituents: first principles density functional theory 
calculations (DFT), as implemented in the QuantumEspresso package [19] and 
an in-house developed Boltzmann transport code [18, 20, 21, 22] to calculate the 
thermoelectric properties. 

First, the band structure of the strained and unstrained Si was calculated using 
the general gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) 
flavor of exchange correlation functional [23] . Fully relativistic and norm-conserving 
pscudo potentials [24] were used to treat the spin-orbit splitting of the Si valence bands 
appropriately. The calculations for the bulk Si were performed with the rhombohedral 
experimental lattice constant ao = 5.434A/V2 for a rhombohedral two atom unit 
cell, which is sketched in the inset of figure 1(c). The strain in [11 Indirection under 
constant volume is simulated by changing the lattice constant a and with that angle 
a r . Throughout the paper the biaxial strain will be given in units of the relative 
change of the in-plane lattice constant, that is the nearest neighbour distance in the 
[lll]-planc as Aa = arm] following the notation of previous works [18, 21, 25, 26]. 

The angle a r is given by cosa r = 1 — ^^(l+Aa/aa) ' That means, tensile in- plane 
strain considers changes Aa /a > and a r > 60°, while compressive in-plane strain 
means Aa /a < and a r < 60°. As used previously in literature, tensile strain along 
[lll]-direction coincides with compressive in-plane strain as denoted here. 

By this trigonal deformation an atomic relaxation of the atomic positions inside 
the unit cell is possible, as the displacement of the two sublattices along [11 Indirection 
is no longer given by symmetry. To obtain the atomic positions of the strained silicon 
we performed structural relaxations using VASP [27]. The atomic positions were 
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optimized such that the Hellmann-Fcynman forces on them were below 0.1 me v/A. 
At the same time the given deformed lattice parameters were not allowed to relax 
and conservation of the unit cell volume was assumed. A volume relaxation at the 
maximum strain Aa /a ±l% resulted in a volume reduction by 0.2% which corresponds 
to a lattice constant change by less than 0.1%. 

As expected, our DFT calculations underestimate the size of the band gap at zero 
temperature and do not reproduce the temperature dependence of the gap. For this 
purpose we implemented a temperature-dependent scissor operator [28], so that the 
strain- and temperature-dependent energy gap E g becomes 



E g (T,— ) = Eg (T = : — ) + U GGA (2) 
a a 



2 



aT 
~ T + /3' 

where E g (T = 0,-^) is the zero temperature gap obtained by our self- 
consistent DFT calculations, C/qga = 0.57 eV is a static correction to reproduce the 
experimental low temperature gap and the third part of Eq. 2 is the correction of 
the temperature dependence of the band gap in a wide temperature range [29], with 
a = 4.73 x 10~ 4 cV /k, T the absolute temperature and (3 = 636 K for bulk silicon. 

Converged results from the first step are basis to obtain the thermoelectric 
transport properties by solving the linearized Boltzmann equation in relaxation time 
approximation (RTA) [20] . Boltzmann transport calculations for thermoelectrics have 
been carried out for quite a long time and show reliable results for metals [30, 31] 
as well as for wide- and narrow-gap semiconductors [18, 21, 32, 33] in the diffusive 
limit of transport. Here the relaxation time is assumed to be constant with respect 
to wave vector k and energy on the scale of ksT. The constant relaxation time 
allows for the calculation of the thermopower S without any free parameter. To 
reproduce experimental findings we parametrized doping dependent relaxation times 
from mobility measurements on unstrained silicon according to ref. [34] by 

r(A) = (( - c/„ . arctan[a • lg( N /N )] + ^2) 

+ 2(lg(^M) 2 ) ■ d, (3) 

with a = 1.8, A = 10 17 cm~ 3 , Ax = 10 17 5 cm- 3 , c = 1500, d = 0.15 fs (a = 1.3, A = 
10 16 8 cm- 3 , Ai = 10 17 5 cm- 3 ,c = 550, d = 0.13 fs) for electron (hole) doping and 
charge carrier concentrations of A between 10 14 cm -3 and 10 22 cm -3 . Nevertheless, 
we state that our relaxation time is not strain-dependent, while it is known, that 
under strain the dominant scattering process alters: For unstrained Si, the room- 
temperature scattering is dominated by optical phonons, i.e., intervalley scattering, 
whereas for strained Si, the scattering by optical phonons is reduced [35, 36]. 

With the transport distribution function (TDF) as termed by Mahan and Sofo 

[37] 



(^3E/^ 3k (<^ii)) 2 ^-^ v 



E=E£ 

(4) 
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the temperature- and doping-dependent electrical conductivity a and thcrmopower S 
are defined as 

™ 1 C l pJu,T) 

a ±Al =e^,T) S ±Al =±^^ y (5) 

Here || denotes the in-plane direction parallel to the a-axis and _L the cross-plane 
direction. The electronic part to the total thermal conductivity accounts to 

^ = t[ C W,T)- 4 o) (M)T) j- (6) 

denotes the band structure of band v, v£ the group velocity and /(^t) the 
FERMi-DiRAC-distribution with chemical potential fi. The chemical potential fi at 
temperature T and extrinsic carrier concentration N is determined by an integration 
over the density of states n(E) 

VBM n+AE 

N = J dEn(E)[f ( ^ T) -l} + J dEn(E)f^ T) , (7) 

(i-A£ CBM 

where CBM is the conduction band minimum and VBM is the valence band maximum. 
The necessary size of AE will be discussed below. 

In a recent work we showed, that the determination of surface integrals in 
anisotropic Brillouin zones is demanding with respect to convergence of the transport 
property anisotropy [22]. So, the constant energy-surface integrations, which are 
required in equation 4, are performed within an extended tetrahedron method 
[38, 39, 40] interpolating the calculated Eigenvalues on an adaptive k-mcsh 
corresponding to a density of at least 44000 k points in the irreducible part of 
the Brillouin zone. £^ll (E, T = 0) was determined on a dense energy mesh with 
a step width of 1 mcV. At vanishing strain Aa /a = the numerical errors of 
o~,,/cr ± and S,./S ± were constantly below 0.1%. In the limit of low carrier concentrations 
N < 1 x 10 14 cm -3 and for larger carrier concentrations in the bipolar conduction 
regimes at high temperatures, convergence of the integrals 4 and 7 was achieved with 
an adaptive integration method for AE of at least lOk^T. 



3. Thermoelectric transport 

3.1. [Ill] -strained silicon: electron doping 

In figure 1(a) the uncorrected band structure of bulk silicon in the rhombohcdral 
unit cell is shown for the unstrained case (black solid line), under 1% tensile strain 
(red dashed line) and for —1% compressive strain (blue dotted line) on relevant 
high symmetry lines. Unstrained silicon has an indirect band gap with conduction 
band minima (CBM) near the F high symmetry point. The CBM consists of sixfold 
degenerate (A 6 ) prolate spheroidal isoenergetic surfaces along six equivalent Y — F 
directions. Due to symmetry of the lattice distortion in [11 Indirection this degeneracy 
holds under applied strain in contrast to strain applied along [001]-direction [18]. 

In the unstrained case, for each Ag-valley the effective masses along the major 
and the minor axis are Mi = 0.91too an d M2 = 0.19mo, respectively. As can be seen 
from figure 1(b) the idea of an effective mass determined by a second order polynomial 
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Figure 1. (color online) Band structure of bulk silicon in the rhombohcdral unit 
cell. In the unstrained case, the F point coincides with the X point of the fee 
Brillouin zone, (a) The bands for the unstrained case (black solid line), under 
1% compressive strain (blue dotted line) and under 1% tensile strain (red dashed 
line) are shown. In (b) a zoom near the conduction band minimum elucidates 
the nonparabolicity of the bands under applied strain, (c) Directional effective 
conduction band masses for silicon under [lll]-strain. The insets show in detail 
the behaviour of the averaged effective mass M on the applied strain, as well as 
the rhombohcdral unit cell [41] (red lines). More details are given in the text. 

fit (dashed lines in figure 1(b)) is valid for the unstrained case, but band warping leads 
to deviations for silicon under strain already for small band occupations. Here, a fifth 
order polynomial fit (solid lines in figure 1(b)) is necessary to reproduce the band 
dispersion, which occurs under applied biaxial [lll]-strain. It is therefore advisable 
to go beyond a simplified effective mass model. Contrary to the conduction bands, 
the [lll]-strain leads to a splitting of the degenerate valence bands, the heavy-hole 
(HH) and light-hole band (LH), similar to the case of [001]-strain [18, 42]. While 
the spin-orbit-split-off band is 40 meV away from the band edge, the HH band lifts 
up energetically under tensile strain, while the LH band lowers in energy. This 
picture reverses under opposite strain conditions [43]. While the indirect gap closes 
linearly under tensile in-plane strain from 0.58 cV in equilibrium to about 0.45 eV at 
Aa /a Q = 1%, the gap size is almost constant within 0.01 eV under compressive strain 
within the considered range [26, 44, 45] (comp. figure 1(b)). The direct gap at T 
decreases slightly under applied strain, more pronounced under compressive in-plane 
strain. 

The influence of biaxial in-plane strain on the fitted effective masses is summarized 
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Figure 2. (color online) Anisotropic thermoelectric transport properties of 
Si for fixed temperature and electron doping concentrations in dependence on 
compressive and tensile strain in [lll]-dircction. Left panels (a)-(c)) correspond 
to an electron doping of 2 X 10 — 8 e/atom (N = 1 X 10 15 cm -3 ) at a temperature 
of 100 K, while right panels (d)-(f) refer to an electron doping of 0.04 e/atom 
(TV = 2 X 10 21 cm~ 3 ) at a temperature of 900 K. On the left axis of each figure 
the relative value compared to the unstrained case is shown, while on the right 
axis the absolute values are given. 

in figure 1(c). The change of the effective mass Mi along the major axis is almost 

symmetric to the applied strain and increases up to 130% of the unstrained value. 

In contrast, the transverse effective masses behave drastically different. While the 

transverse effective mass M2 increases under tensile strain, the effective mass M3 

perpendicular to M2 decreases. This behaviour is reversed under applied compressive 

strain. These results are in good agreement to other findings [44, 45], but add up 

to GW calculations were only one transverse mass with an almost constant value 

3 

was found [26]. The strain-dependent averaged effective mass M = n (M^ 1 / 3 , often 

i=l 

referred to as density-of-states effective mass, is shown additionally as an inset in 
figure 1(c). 

In the following the influence of trigonal distortion on the thermoelectric transport 
of Si under electron doping will be discussed. For this purpose two doping and 
temperature regimes are considered. The first, at a low temperature of T — 100 K 
and low charge carrier concentration of N = 1 x 10 15 cm -3 , is suitable for metal- 
oxide-semiconductor device applications. The results for the electrical conductivity, 
thermopower and power factor in dependence on the in-plane strain are shown in 
figure 2(a), (b) and (c), respectively. Under tensile strain the in-plane electrical 
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conductivity increases almost linearly, while the cross-plane component a ± decreases 
almost comparable. For compressive strain the behaviour reverses, with the cross- 
plane conductivity being enhanced up to 23 % at 1% compressive strain, while the 
in-plane conductivity decreases to nearly 23% of the unstrained value under 1% tensile 
strain. In the limit of a degenerate semiconductor at low temperatures and small 
charge carrier concentrations this results can be completely understood within an 
effective mass calculation [22] . With noticeable variation of the electrical conductivity 
under applied strain, the thcrmopower is almost unaffected. At low temperatures 
only a small energy window near the band edges plays an important role for the 
determination of the thermopowcr. As the functional change of the coefficients 
(/•*' * s determined by parabolic bands and a strain-dependent shift of the 
chemical potential, a strong change in the thermopower can not be expected. However, 
a slight upward tendency for the thermopowcr under compressive and tensile strain can 
be stated. As S oc (-Ecbm/vbm — m) oc M, this can be directly linked to the strain- 
dependence of M (cf. figure 1(c)). For tensile strain an anisotropy of the thermopower 
is apparent, which seems to be suppressed for compressive strain. This could be linked 
to stronger deviations from the isotropic effective mass for tensile strain as shown in 
figure 1(b). 

The strain-dependence of the resulting power factor PF is shown in figure 2(c). 
Due to the weak impact of the thermopower, the behaviour of the power factor is 
dominated by the electrical conductivity dependence on the applied strain. At low 
temperatures and small charge carrier concentrations, such a behaviour has already 
been observed for biaxially strained silicon [18]. However, the additional power-output 
described by the power factor is enhanced by 30% under applied 1% compressive strain 
for PF^and by 25% under applied 1% tensile strain for PF |r We note, that this low 
temperature and small doping case is not feasible for thermoelectric power generation, 
but could give inside in parasitic effects which play a role in metal-oxidc-semiconductor 
devices. 

Usual conditions for silicon-based thermoelectric applications, such as high 
temperature of 900 K and large charge carrier concentrations N = 2 x 10 21 cm -3 are 
assumed in figures 2(d)-(f). At temperatures of 900 K the electronic band structure 
on a width ±AE = ±800 meV around the position of the chemical potential has 
to be considered, which makes a simplified description of the electronic transport 
properties within a spherical band picture inaccurate. However, the dependence of the 
electrical conductivity (ref. figure 2(d)) on the applied [lll]-strain is almost preserved, 
even if accidentally. This is in contrast to biaxial strain in [001]-direction, were the 
strain-dependence induced by reoccupation of bands was suppressed under higher 
temperatures and dopings [18, 35]. 

In figure 2(e) the anisotropic thermopower under trigonal distortion is shown. 
The cross-plane thermopowcr shows a monotonous increase from 85% to 110% of the 
unstrained thcrmopower value of S ± = 115/xV/K varying from compressive to tensile 
strain. The in-plane component S„ decreases to 85% of the unstrained thermopower 
value under 1% tensile strain. A very similar behaviour was found for the thermopower 
of biaxially strained silicon in [001]-direction [18]. The compensation effects of the 
enhanced thermopower and decreased electrical conductivity and vice versa, are well 
known for thermoelectrics under strain [18, 21, 46]. This scenario holds for [111]- 
straincd silicon, too. In figure 2(f) the shown anisotropic power factor in in-plane 
and cross-plane direction is always smaller than the power factor of the undistorted 
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system, while PF^ being at least stable under small values of compressive strain. 
We mention that the absolute values (cf right scales in Figs. 2(a) and (d)) of the 
electrical conductivity are increased remarkably compared to the low doping case as 
expected. Due to this, the absolute value of the power factor raises in its absolute 
value, but unfortunately does not show an enhancement due to mechanical strain in 
[lll]-direction. Furthermore, compared to the low-doping/low-tcmpcrature regime 
(figure 2(c)) the power factor does not show noticeable anisotropy between in-plane 
and cross-plane components. 

3.2. [Ill] -strained silicon: hole doping 

As it is well known, thermoelectric devices use two types of semiconductors, namely 
n-type and p-type, which are connected in series [6]. Therefore, the influence of 
biaxial [lll]-strain on hole-doped silicon is presented in figure 3 in the same way 
as done for the electron-doped case. In the low-doping/low-tcmperature regime an 
enormous enhancement for the cross-plane electrical conductivity asunder sufficient 
tensile strain can be found (cf figure 3(a)), while the in-plane component cr.. decreases 
more slightly under the same strain conditions. This behaviour can be linked to a 
changed subband structure. As mentioned before, strain lifts the degeneracy of LH 
and HH bands around the T point and alters the curvature, that is the effective mass, 
of both bands. Under applied strain, the valence bands become highly anisotropic and 
a crossover between bands occurs so that they even loose their original LH and HH 
meaning[18, 25]. Extended discussions here on can be found in refs. [43, 47]. 

The thcrmopowcr of p-type silicon is shown in figure 3(b). The thermopower 
in-plane and cross-plane decreases slightly under tensile, as well as under compressive 
strain. The anisotropy of the thermopower is moderate. As previously reported 
[18, 48], the thermopower depends strongly on the number of occupied carrier pockets. 
A higher valley degeneracy at a fixed charge carrier concentration leads to an increased 
thcrmopowcr. As already mentioned the formerly degenerate HH and LH bands split 
under tensile and compressive strain. At low hole concentrations and low temperatures 
only the former HH (LH) band is occupied. This leads directly to a reduction of and 
S ± . Consequently, the accompanied power factor (cf. figure 3(c)) is also reduced in 
its maximal possible enhancement, but follows in principle the behaviour given by the 
electrical conductivity. 

For high temperatures and high hole concentrations the results are shown in 
figures 3(d)-(f). As the carrier concentration is raised by nearly six orders of magnitude 
with respect to the low doping case, the absolute value of the electrical conductivity 
raises in the same order (note the right scale of figure 3(d)). Unfortunately, due to the 
higher band occupation and the broader smearing of the Fermi-Dirac- distribution, 
effects of redistribution in strain-split bands do not play a role any more. Moreover, the 
effects of reduced and increased effective masses cancel each other leading to negligible 
change, in absolute values as well as in anisotropy, of the hole electrical conductivity. 

For the hole thcrmopowcr shown in figure 3(e) this behaviour is still valid. For 
the in-plane thcrmopowcr S,, no significant influence of either compressive or tensile 
strained could be found. A minor dependence on applied [lll]-strain is observed 
for S ± . Here, the cross-plane thermopower decreases from 100 //V/K to 91/kV/K for 
strain values varying from 1% compressive to 1% tensile strain. 

Comprising the results for the electrical conductivity and thermopower, the 
resulting power factor under hole doping is shown in figure 3(f). No evident influence 
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Figure 3. (color online) Anisotropic thermoelectric transport properties of Si for 
fixed temperature and hole doping concentrations in dependence on compressive 
and tensile strain in [111 ]-direction. Left panels (a)-(c)) correspond to a hole 
doping of 2 X 10" 8 e/atom (N = 1 X 10 15 cm" 3 ) at a temperature of 100 K, while 
right panels (d)-(f) refer to a hole doping of 0.04e/atom (N = 2 X 10 21 cm -3 ) at a 
temperature of 900 K. On the left axis of each figure the relative value compared 
to the unstrained case is shown, while on the right axis the absolute values are 
given. 



of [lll]-strain on the power factor could be found in the thermoelectrically relevant 
temperature and hole doping regime. Furthermore, the absolute value of the power 
factor is about 3-4 times smaller than in the comparable electron doped case, which 
is mainly caused by the higher hole scattering rate as assumed in Eq. 3. At varying 
strain not only the absolute values of the thermoelectric properties change, but also the 
optimal charge carrier concentrations to obtain these maximized values. For biaxially 
[001]-strained silicon it was found, that the optimal doping range can change by a 
factor of two, while reducing the power factor up to 20% if not adapting the charge 
carrier concentration [18]. Therefore, the in- plane and cross-plane power factor under 
varying electron- and hole-doping, as well as varying [lll]-strain is shown in figure 4. 
The maximized power factor at optimal charge carrier concentration at a given strain 
state is emphasized by a black-dotted line in figures 4(a)- (d). Comparing silicon 
under electron doping (cf figures 4(a), (b)) and hole doping (cf. figure 4(c), (d)), it is 
obvious, that electron-doped silicon shows a much stronger variation of the optimal 
doping range. For the in-plane component PFnand electron-doping the optimal 
carrier concentration decreases by about a factor 5, from N = 5 x 10 21 cm -3 to 
N = 1 x 10 21 cm -3 with the strain changing from 1% compressive to 1% tensile strain. 
For the corresponding cross-plane power factor PF^the optimal carrier concentration 
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Figure 4. (color online) In-plane and cross-plane power factor of Si at fixed 
temperature of 900K under electron (a,b) and hole doping (c,d) for varying charge 
carrier concentration and applied [lll]-strain. The black circles emphasize the 
position of optimal doping at a certain strain state to maximize the power factor. 
Note the different scale for electron and hole doping. 



is about N = 1 X 10 cm~ 3 at 1% compressive strain and increases by about a factor 5 
under same strain conditions. At hole-doping no evident change of the optimal doping 
could be found at varying strain, while the absolute value depends only weakly on the 
applied strain (cf. figure 3(f)). 

To summarize our findings, sections along the path of optimal electron doping 
are shown in figure 5. With this, it is obvious that even with optimized doping 
no enhancement by tensile or compressive [lll]-strain can be obtained for the in- 
plane power factor PF|| . For the cross-plane component PF^ , which is more relevant 
for possible SL's, an increase of the power factor of about 4% at 1% compressive 
strain was found. To obtain this rather small enhancement the electron charge carrier 
concentration has to be reduced by about a factor of 2.5 compared to the unstrained 
case. 



3.3. Strained Si/Ge-SL on Si[lll] 

By introducing the concepts of carrier pocket engineering [15, 16, 17] and phonon- 
glass/electron-crystal [49, 50] to semiconducting SLs an enormous leap forward to 
maximize the thermoelectric FOM was proposed. Indeed, several proofs-of-principlc 
showed a remarkable enhancement of the FOM for thermoelectric semiconducting 
heterostructures [7, 8, 17, 51]. With the thermal conductivity of SL's far below 
their alloy limit [52, 53, 54, 55] and their constituents bulk values, a main task 
in optimizing the FOM is to enhance or at least to retain advantageous electronic 
properties of the bulk materials, so to say the power factor. For silicon-based 
SL's, carrier pocket engineering can be triggered by lattice strain. Using Sii^/Ge^ 
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Figure 5. (color online) Maximized in-plane and cross-plane power factor of 
[lll]-straincd Si at fixed temperature of 900K under optimal electron doping. The 
figure highlights the path of the black circles in fig. 4(a) and (b). N = 0.01 e/atom 
corresponds to N = 5 X 10 20 cm -3 . 



substrates in [111] and [100] orientation a ZT = 0.96 and 0.24, respectively, were 
predicted for strain-symmetrized§ Si(20A)Ge(20A) SL's at room temperature. The 
latter case was experimentally confirmed with ZT = 0.1 at N as 1 x 10 19 cm~ 3 
[17], which is nevertheless about a sevenfold enhancement relative to bulk Si 
[13, 18]. These experimental findings encourage further research for strain-non- 
symmetrized Si(20A)/Ge(20A) SLs in [001]-orientation and Si(15A)/Ge(40A) SLs in 
[lll]-orientation, with ZT = 0.78 and 1.25, respectively, predicted at T = 300 K [15]. 

While being in principle possible for very thin films [58, 59], to our best knowledge 
state-of-the-art thin film technology docs not enable strain-non-symmetrized SL's with 
satisfactory structural qualities and thicknesses for thermoelectric applications so far 
[56, 57, 60, 61]. Nevertheless, in the following the thermoelectric transport properties 
of a strain-non-symmetrized Si(5A)/Ge(5A) SL in [lll]-orientation will be discussed, 
as the largest enhancement of the FOM is expected here[15]. 

The used Si(5A)/Ge(5A) SL is represented by a hexagonal six-atom unit cell 
(ref. inset in figure 7(a)) with point group symmetry C% v and a fixed in-plane lattice 
constant of a = 5.434A was used to simulate the bulk silicon substrate. Structural 
optimization of the atomic positions and the c-axis elongation of the unit cell were 
obtained using VASP [27]. The distinct interlayer distances 6 in [11 Indirection 
were determined as 5si 1 -si 2 — 2.359A, Ssi 2 si 3 — 0.784A, Ss^-Ga = 2.409A , 
Sce 4 -Ge 5 = 0.929A, SGe 5 ~Ge 6 = 2.470A, (fcee-Su = 0.85lA and are in good agreement 
with previous calculations [62]. The ratio c /a = 2.551 shows an increase of the lattice 
constant in c-dircction by about 4% compared to bulk Si. This is clearly dictated by 
the fixed in-plane Si lattice constant and the larger volume of Ge. Due to this, the Ge 
layer can be seen as compressively strained in [11 Indirection. 

As has already been mentioned in literature [62, 63], the face centred cubic (fee) L 
high symmetry point in the [111] direction folds onto the hexagonal A high symmetry 
point at a c /a-ratio of 2.449. In addition, the fee X point is equivalent to the hexagonal 
M point, while the symmetry directions fee FX and hexagonal FM are inequivalcnt. 
This is due the fee X point lying in an adjoining Brillouin zone. Under biaxial [111]- 
strain two inequivalcnt sets of the eight L points occur. There are two points along 
c-axis in [111] direction of growth which fold onto the A point and those in the six 

§ For an introduction to strain-symmetrized and strain-non-symmetrized Si/Ge supcrlattices we refer 
to the publications of Kasper et al. [56, 57]. 
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Figure 6. (color online) (a) Band structure of the Si/Ge SL. The color code of the 
bands refers to the atomic character of the bands. Red dots refer to wavefunction 
character of pure germanium, while blue dots refer to pure silicon like character, 
(b) band structure of silicon (blue, dashed lines), germanium (red, dotted lines) 
and Si/Ge (black, solid line). All three configurations are calculated in the lattice 
of Si/Ge. To allow comparison with (a), the bands are energetically matched at 
the CBM. The figure has not been corrected for the band gap error. 



directions equivalent to [111] hereafter noted as L. 

In figure 6 the band structure for the fully relaxed Si/Ge-SL in the hexagonal 
unit-cell is shown. While in figure 6(a) the site-resolved bands for the Si/Ge-SL arc 
depicted, in figure 6(b) the bands for pure Si on the fully- relaxed Si/Ge-positions are 
shown (blue, dashed lines), as well as all sites occupied by Ge (red, dotted lines) and 
the bands referring to the original Si/Ge-SL bands as shown in figure 6(a) (black, 
solid lines). Two main insights can be drawn. First, the VBM is located around the T 
point and has an almost pure Ge character, given by the red dots in figure 6(a). This 
is mainly due to the smaller band gap of Ge and the accompanied band offset between 
Ge and Si. Furthermore, the compressive biaxial strain in the Ge layers decreases the 
gap size and favours a direct band gap at T instead of an indirect one between T and 
L [64. 45, 65]. Valence bands of mixed Si/Ge character come in to play around 0.65 eV 
below CBM and suggest the local indirect Si band gap between V and M to be almost 
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retained bulk like. Second, the CBM is located at the M point and shows a strong 
mixing of Si and Ge character. Furthermore a strain-induced lifting of degeneracy 
occurs at the M point lowering a band of mixed character to the band edge and lifting 
a Si-like band upwards in energy. As only the Ge layers in the SL are compressively 
strained while the Si layers are nearly unstrained, the Ge L-point CBM valleys split 
into valleys located at A higher in energy and in L valleys lower in energy [15, 65, 66]. 
These findings on the CBM and VBM characteristics are in agreement to experimental 
studies [59]. With that the uncorrected band gap decreases to around 178 meV, which 
is about 30% of the uncorrected GGA gap for unstrained bulk Si. Furthermore, the 
effective masses at the CBM lower to Mi=0.125mo, M 2 =0.026m and M 3 =0.010m 
|| . Applying an effective mass approach [22] we find the conductivity anisotropy for 
energies near the CBM to <T ||/o i =0.6, which clearly prefers cross-plane transport under 
electron-doping. 

The energy-dependent transport distribution function in the in-plane and cross- 
plane direction and their ratio are displayed in figure 7(b) for the Si/Ge-SL and 
isotropic unstrained bulk Si. For the same systems the density of states are shown in 
figure 7(a). As can be seen from figure 7(b) unfortunately the conductivity anisotropy 
near the valence band edge strongly increases to values of 8 around 0.58 eV below 
CBM, clearly suppressing cross-plane electronic transport under p-typc doping. This 
behaviour is largely due to the localization of the Ge like VBM states, in space as well 
as energy. As can be deduced from figure 6(a) bands in the cross-plane direction (here 
TA) show pure Ge band character in an energy range of 0.18 — 0.74 eV below CBM. 
Clearly, in this energy range cross-plane conduction is suppressed by states localized 
in the Ge layers and almost vanishing in the Si layers leading to a strongly increased 
conductivity anisotropy. At an energy of 0.38 eV below CBM at T a light band with 
strong Si/Ge mixed character appears (cf. black solid line in figure 6(b)) leading to 
cross-plane transport through the Si and Ge layer forcing <J n/<y ± to decrease, while 
saturating for values above 1, still indicating a preferred in-plane transport under 
hole-doping. 

In figure 8(a),(b) the doping-dependent thermopower and power factor for the 
Si/Ge-SL are shown, respectively. As a comparison the reference values of bulk silicon 
are stated as black, dashed-dotted lines. We note, that a temperature dependency of 
the energy gap was introduced applying Eq. 2. For the Si/Ge-SL the parameters 
are chosen as Ugga = 0.78 eV, a = 4.76 x 10 _4eV 7 K and f3 = 395 K. However, in the 
highly degenerate limit (N > 1 x 10 20 cm" 3 ) the temperature dependence of the gap 
plays a negligible role even for temperatures above 900 K. 

From figure 8(a) it can be seen, that the thermopower in the Si/Ge-SL under 
electron doping (blue lines in the lower panel) is comparable to bulk silicon and follows 
a Pisarenko relation [68]. Under hole doping (red lines in the upper panel) the 
thermopower is suppressed compared to bulk silicon by about 80 /iV/K for S n and 
S ± in the relevant doping regime, which might be linked to a changed functional 
behaviour of the TDF cf'V(iJ,,T). The latter can be deduced from figure 7(b) where 
apparent differences in the functional behaviour of cr ± .. , which is proportional to 
^±\\ ■^')' arc visible especially in the valence bands of bulk Si and the Si/Ge-SL. The 
clear deviation of S^from the Pisarenko relation in the vicinity of N as 3 x 10 21 cm" 3 

| The referring Eigenvectors were determined as e\ = (0.8,0,0.6), e2 = (0,1,0) and e% = 
(0.6,0, -0.8). 

If Here we adapted experimental data for Sio.sGeo.s alloys from Ref. [67]. 
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Energy E-E c (eV) 

Figure 7. (color online) (a) Density of states for bulk silicon (green, dashcd- 
dotted line) and the Si/Ge SL (black, solid line). As an inset the hexagonal 
unit cell of the Si/Ge SL is shown [41]. Furthermore the interlayer distances are 
labelled, as referred to in the text, (b) Electrical conductivity in dependence on 
the position of the chemical potential fi at zero temperature, shown for bulk silicon 
(green, dashed-dotted line) and the Si/Ge SL in the in- plane (red, solid line) and 
cross-plane direction (blue, solid line). The conductivity anisotropy (black, dashed 
line referring to the right axis) is stated for the Si/Ge SL. The cross at the CBM 
is the value obtained from an analytical effective mass approach. 

is related to the strongly increased conductivity anisotropy 0.58 eV below VBM 
as shown in figure 7(b). Here, the strong suppression of <j ± causes a larger S ± . 
Nevertheless, this slight enhancement of S ± is not reflected in the power factor of the 
p-type Si/Ge-SL. As shown in figure 8(b) S,, , as well as S^are always smaller than the 
values for bulk Si under hole doping (red lines). Obviously, the suppressed electrical 
conductivity, especially in the cross-plane direction, is responsible for this result. A 
power factor of about 9 /iW/cmK 2 is found for PF,, at N ps 3 x 10 20 cm~ 3 , while in the 
cross-plane direction the same value can be stated at huge values of ./V ss 3x 10 21 cm -3 , 
clearly evoked by the anomaly in the thermopower. More interesting is the case of 
electron doping (blue lines in figure 8(a)). With the thermopower's behaviour almost 
bulk like and conductivity anisotropies a u/<r ± below 1 an enhanced power factor in 
the required cross-plane direction is found. Compared to bulk silicon the PF^is 
enhanced by 10% and reaches a value of 60 /iW/cmK 2 at an electron concentration of 
N ps 7 x 10 20 cm~ 3 . With that the optimal charge carrier concentration is four-times 
smaller compared to bulk Si. For the in-plane component PFn almost no reduction can 
be seen, while the maximal value of 53/xW/cmK 2 is shifted to slightly larger charge 
carrier concentrations. 
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Figure 8. (color online) (a) In-plane (solid lines) and cross-plane (dashed lines) 
doping-dcpendent thermopowcr at 900 K for the Si/Gc SL under electron (thick, 
blue lines) and hole-doping (thin, red lines). For comparison the values for bulk 
silicon are given (black, dashcd-dottcd lines), (b) In-plane (solid lines) and cross- 
plane (dashed lines) doping-dcpendent power factor at 900 K for the Si/Ge SL 
under electron (thick, blue lines) and hole-doping (thin, red lines). Again, for 
comparison the values for bulk silicon are given (black, dashcd-dottcd lines). 
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Figure 9. (color online) (a) Doping-dependent cross-plane figure of merit 
ZTx of the Si/Ge SL (blue lines) and bulk Si (black lines) under electron- 
doping at different temperatures. (b) the total thermal conductivity k± in 
cross-plane direction for the Si/Ge SL (blue lines) and bulk Si (black lines) 
under electron-doping at different temperatures. While the electronic part ft e !X 
was calculated, the lattice part K p hx was estimated from experiments [13, 55] 
and is constant for varying charge carrier concentration. (c) the calculated 
cross-plane Lorenz function Lx = K el± ■ (ffiT)" 1 related to the metallic limit 
L = 2.44 X lO" 8 WO/K 2 . 
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3.4- Towards the figure of merit 

In figure 9(a) the FOM in cross-plane direction for the Si/Ge-SL (blue lines) and the 
more promising electron-doped case is shown. Different temperatures are chosen to 
demonstrate the evaluation of maximal ZT and the range of optimal charge carrier 
concentration. As a comparison the FOM for bulk Si is shown as black dashed-dotted 
line, too. To present results for the FOM, knowledge of the thermal conductivity is 
relevant. For this purpose, the electronic part of the thermal conductivity K e i± was 
calculated applying Eq. 6, while the lattice part n p h±_ was taken from experiment. 
Here, n p h± = 36W/mK [13] was used for bulk Si at 900 K, while k p h± — 1.9W/mK 
was used for the Si/Ge-SL [55]. With the latter being rather optimistic, since smaller 
than the expected nano-alloy limit of 2.5W/mK [69], but achievable in Si/Ge-SL 
[53, 55, 70]. However, from figure 9(b) one can deduce that for thermoelectric 
reliable charge carrier concentrations above 4 x 10 20 cm -3 the electronic contribution 
to the total thermal conductivity dominates over its lattice part. At 900 K for 
N ~ 1.5 x 10 20 cm -3 the electronic part amounts to 0.73 W/mK clearly smaller than 
the lattice part. This contribution increases significantly at higher temperatures and 
charge carrier concentrations omitting higher absolute values of the FOM. However, 
ZTj_ above unity can be reached for operating temperatures higher than 750 K and 
electron-doping N ss 1.5 — 3 x 10 20 cm -3 . 

At room-temperature a ZT± ps 0.2 is achieved in a broad doping range of 
JV »i 5 — 25 x 10 19 cm -3 . Besides this value being an order of magnitude higher 
than the bulk Si value of ZTw 0.01 [13, 15, 18], it is still less than the postulated 
values of ZTj_ = 0.96 and ZTj_ = 1-25 by Koga et al.[15] for strain-symmetrized 
and strain-non-symmetrized Si/Ge-SL, respectively. With that, and using a rather 
conservative value of K p hi_ = 7.3 W/mK for their estimations, they expect enormous 
power factors of PF^ w 250^W/cmK 2 and PFj^ rj 340/iW/cmK 2 for the strain- 
optimized Si/Ge-SL in [lll]-direction. We found PF ± « 15^W/cmK 2 at 300 K and 
N ps 1 x 10 20 cm" 3 for the Si/Ge-SL. Even though a convergence of carrier pockets is 
not fully achieved in our superlattice, power factors beyond 200 ^W/cmK 2 seem to be 
very sophisticated, as state-of-the art power factors near or above room-temperature 
are well below 100 /iW/cmK 2 [7, 71, 72]. A benefit from thermoionic emission even at 
moderate temperatures can not be expected in Si/Ge-SL [73, 74, 75]. We note, that 
at very low temperatures below 10 K huge PF of about 100 — 1000/iW/cmK 2 were 
reported for bulk Fe 2 Sb 2 and ZnOi_ K Sc K [76, 77]. 

As an addition in figure 9(c) the doping-dependent Lorenz function Lj_ = 
K e i± ■ (cr±T) _1 as defined via equations (5) and (6) is presented. From figure 9(c) 
it is obvious that the Lorenz number L± can be substantially different from the 
metallic limit L$. Nevertheless, for very large charge carrier concentrations and 
the chemical potential located deep inside the conduction band, L± almost coincides 
with Lq. At intermediate and thermoelectric relevant charge carrier concentrations 
of A rj 5 — 50 x 10 19 cm -3 L± can be much smaller than Lq. For 900 K and 
N w 1.25 x 10 20 cm~ 3 we find a minimal value of L± w 0.685Lo. At smaller charge 
carrier concentrations L± rapidly increases and reaches L± ss 7Lq for electron charge 
carrier concentrations of A < 1 x 10 18 cm -3 in the intrinsic doping regime. At 
decreasing temperatures minimal values of the Lorenz function are obtained at much 
smaller charge carrier concentrations. Furthermore, the maximal values of L shifts 
to smaller charge carrier concentrations, too, and can reach huge values of L at very 
low temperatures and charge carrier concentrations. The effect, which is responsible 



19 



for the suppression of the Lorenz function to values below the metallic limit L is 
termed bipolar thermodiffusion effect [78, 79, 80] and maximizes for positions of the 
chemical potential near the band edges. However, a Lorenz function L =/= Lq can 
have consequences for the determination of the thermal conductivity. The Lorenz 
factor is generally used to separate n e i and n p h- At thermoelectric advisable charge 
carrier concentrations applying the metallic value Lo to determine the lattice thermal 
conductivity could lead to an ovcrestimation of the electronic thermal conductivity, 
and consequently to an underestimation of the lattice contribution. 

4. Conclusions and Outlook 

With the presented results we have shown, that strain in [11 Indirection is not 
sufficient to significantly enhance the thermoelectric transport properties in bulk Si 
for energy harvesting applications. In the low-temperature and low-doping case large 
enhancements were found at tensile strain for PF|| (electron-doping) and PF^ (hole- 
doping) and under compressive strain for PF ± (electron-doping). This could have 
a negative impact for metal-oxide-semiconductor devices involving [lll]-strained Si. 
Here, in the low-temperature and low-doping regime, small temperature gradients in 
the devices could lead to an additional parasitic electrical power, which can be way 
larger than expected from unstrained bulk. 

Enhancements found in the high-temperature and high doping regime were 
distinct smaller. Here slight enhancements of 5% for PF^were found under 
compressive strain. The more interesting is, that the power factor is robust against 
[lll]-strain, especially under hole-doping. Thus, thermoelectric SLs based on [111]- 
strained Si could provide an enhanced FOM, as n p h is most likely reduced in SLs. 
We note, that due the to high bulk thermal conductivity and the only modest gain in 
the power factor by [111] -strain engineering, bulk silicon remains an unfavourable 
thermoelectric, even if the electronic transport properties are strain optimized. 
However, from comparison with earlier studies on biaxially [001]-strained silicon [18], 
we confirm that strain in [11 Indirection, e.g. in silicon-based SLs, should be preferred, 
as the carrier pocket degeneracy is retained and therefore the thermopower and power 
factor can be maximized. 

To deal with that, we investigated the anisotropic thermoelectric transport of an 
[111] -oriented Si/Ge superlatticc. At a first glance we have shown that no degradation 
of the electronic transport by the heterostructure is expected for electron doping, 
while even showing an enhancement of 10% in PF^ compared to bulk Si. Assuming a 
decrease in lattice thermal conductivity a large enhancement in ZT to 0.2 and 1.4 is 
achieved at 300 K and 900 K, respectively. Under hole doping the electronic transport 
in the Si/Ge-SL is heavily suppressed due to quantum well effects. Here the cross- 
plane power factor PF^is expected to show only around 50% of the bulk maximal 
value, leading to small ZT values. 
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